Mean-field limit of systems with multiplicative noise 
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A detailed study of the mean-field solution of Langevin equations with multiplicative noise is 
presented. Three different regimes depending on noise-intensity (weak, intermediate, and strong- 
noise) are identified by performing a self-consistent calculation on a fully connected lattice. The most 
interesting, strong-noise, regime is shown to be intrinsically unstable with respect to the inclusion 
of fluctuations, as a Ginzburg criterion shows. On the other hand, the self-consistent approach is 
shown to be valid only in the thermodynamic limit, while for finite systems the critical behavior is 
found to be different. In this last case, the self-consistent field itself is broadly distributed rather 
than taking a well defined mean value; its fluctuations, described by an effective zero-dimensional 
multiplicative noise equation, govern the critical properties. These findings are obtained analytically 
for a fully connected graph, and verified numerically both on fully connected graphs and on random 
regular networks. The results presented here shed some doubt on what is the validity and meaning 
of a standard mean-field approach in systems with multiplicative noise in finite dimensions, where 
each site does not see an infinite number of neighbors, but a finite one. The implications of all this 
on the existence of a finite upper critical dimension for multiplicative noise and Kardar-Parisi-Zhang 
problems are briefly discussed. 

PACS numbers: 02.50.Ey,05.50.-a,64.60.-i 



I. INTRODUCTION 

Problems susceptible to be mathematically represented 
by stochastic (Langevin) equations including a multi- 
plicative noise abound not only in physics, but also in 
biology, ecology, economy, or social sciences. In a broad 
sense a Langevin equation is said to be multiplicative if 
the noise amplitude depends on the state variable/s it- 
self/themselves pj. In this sense, problems exhibiting 
absorbing states, i.e. fluctuation-less states in which the 
system can be trapped, are described by equations whose 
noise amplitude is proportional to the square-root of the 
(space and time dependent) activity density, vanishing 
at the absorbing state 0. Systems within this class 
are countless: propagating epidemics, autocatalytic re- 
actions, reaction-diffusion problems, self-organized criti- 
cality, pinning of flux lines in superconductors, etc. 0. 

In a more restrictive sense, the one we will use from 
now on, it is customary to restrict the term multiplica- 
tive noise (MNVto noise amplitudes linear in the activity 
density 0, 0, • Such equations appear ubiquitously in 
economics, optics, population dynamics, study of insta- 
bilities, etc. In all these cases, multiplicative noise terms 
appear rather straightforwardly when constructing (more 
or less rigorously) stochastic representative equations. In 
their spatially extended version, Langevin equations in- 
cluding a MN (sec Eq. below) were first proposed, to 
the best of our knowledge, in the context of synchroniza- 
tion of coupled map lattices 0, and soon after studied 
in Refs. |1 H E3, [H El • Such equations describe dif- 
ferent situations; in particular, there has been a recent 



interest in their application to non-equilibrium wetting 
|13| and also to synchronization problems in extended 
systems ^Hi see Ref- H3 f° r a recent review. In all 
these cases, there is a phase transition between an active 
phase in which the system has a stationary non- vanishing 
activity and another one in which the density field falls 
continuously toward zero without ever reaching it in fi- 
nite time. The analogies and differences between this 
family of models and the one with a square-root type of 
noise, in which the absorbing state is reachable within a 
finite time, have been discussed in [l6lll7| . 

It is remarkable that by using the so-called Cole-Hopf 
transformation, the spatially extended MN Langevin 
equation can be mapped into a non-equilibrium inter- 
face, as represented by the Kardar-Parisi-Zhang (KPZ) 
equation 0, in the presence of a bounding wall ^3 . 
In this language, the active phase describes interfaces 
pinned by the wall while the absorbing one corresponds 
to depinned KPZ-likc interfaces moving away from the 
wall. 

In recent years many aspects of systems with MN have 
been elucidated. For example, a renormalization group 
approach has been constructed, scaling relations derived, 
and the large- TV limit (where N is the number of compo- 
nents) studied It is well established that above two- 
dimensions there are two different regimes; a weak-noise 
and a strong-noise one, in full analogy with the known 
phenomenology of KPZ interfaces [l9|. There is however 
a crucial point which remains to be fully understood: the 
mean-field behavior of such systems. Another interesting 
mapping is that in the absence of the non-linear term, the 
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MN equation corresponds to the equation governing the 
evolution of directed polymers in random media (See [l9l | 
and references therein, as well as [TEjl. 

In Langevin equations with additive noise, as for in- 
stance the Model A describing the universality class of 
Ising like transitions [2(j, the mean- field approximation 
can be obtained in a number of equivalent ways, all of 
them leading to the same results with different degrees 
of sophistication. For example, mean-field critical expo- 
nents can be obtained (i) by discarding the noise and 
solving the remaining deterministic equation, (ii) as the 
most probable solution in a path integral formulation, 
(iii) self-consistently by assuming that each site sees the 
average of the remaining, (iv) by naive power counting 
in the corresponding action, (v) as the lowest order in a 
perturbative loop expansion, etc. 

Contrarily, in systems with MN the meaning of the 
mean-field solution is much less clear. Indeed, early stud- 
ies led to conflicting results depending on the consid- 
ered approximation method. For example, by remov- 
ing the noise, much of the physics is lost and trivial re- 
sults (identical to those for additive noise) are obtained. 
Also, if all spatial dependence is eliminated (by remov- 
ing the Laplacian term in the MN equation), one is left 
with a solvable zero-dimensional equation which does not 
reproduce faithfully the rich mean-field phenomenology. 
Therefore, contrarily to more standard problems, both 
space-dependence and noise have to be retained in order 
to construct a sound mean-field solution. 

If the bounding wall is eliminated then one is left with 
the directed polymer in random media equation, for which 
many results are available. In particular, Derrida and 
collaborators worked out a solution on the Cayley tree, 
mean- field results, 1/d expansions, and solutions on hi- 
erarchical lattices exits |2l|. Also, Mezard and Parisi 
derived a variational approach in replica space p^. and 
also Fisher and collaborators reached also interesting re- 
sults on these issues [23|. For a more detailed review 
on results for directed-polymers see [l9l | and references 
therein. 

On the other hand, restoring the bounding wall (which 
is the case we are interested in) Birner et al. [IJl (see 
also H [13 El) performed a self-consistent calculation 
and reported on the existence of two different mean-field 
behaviors: a weak-noise regime (in this particular case, 
the noise can be completely disregarded) and a strong- 
noise one exhibiting non-universal exponents depending 
on the noise amplitude. This is in agreement with the 
field theoretic expectation of two different behaviors for 
MN-like and KPZ-like equations in high dimensional sys- 
tems (where mean-field results are expected to be valid) 
[HEU. These results are difficult to compare with the 
abovementioned ones obtained for directed polymers in 
random media, which correspond to the deppined phase 
of the full problem, while the self-consistent approach is 
intrinsically devised for the pinned phase. 

Let us also underline that a full understanding in terms 
of path integrals and extremal paths is still missing de- 



spite some efforts in this direction 0, HH ■ 

Together with the determination of the right mean- 
field theory, another relevant and highly debated issue is 
to establish the upper critical dimension d c , above which 
mean-field results hold. As some exponents, as the dy- 
namical one z have been claimed H, to coincide for 
MN and KPZ, both problems are expected to have the 
same upper critical dimension. While there is consensus 
that the mean-field weak-noise regime should be valid 
above d c — 2 (coinciding with the critical dimension for 
weak- noise Edwards- Wilkinson interfaces jljj) there ex- 
ist highly conflicting results (some pointing to d c = 4 j2|J 
and some supporting d c = oo |27jp for the strong- noise 
one regime. 

In any generic problem, for sufficiently high dimen- 
sions every site in a spatially extended system "sees" the 
average of its neighbors (assuming the space has been 
discretized) which can be correctly approximated by its 
mean value (the distribution of the average values is well 
described by the standard central limit theorem above 
the upper critical dimension) defining in this way a sound 
mean-field solution. In the case of MN, as we will show, 
the situation is somehow anomalous: in the strong- noise 
regime the mean-field solution itself breaks down in the 
neighborhood of the critical point once fluctuations are 
taken into account (as a version of the Ginzburg criterion 
shows). This stems from the primary role played by fluc- 
tuations and rare events in multiplicative processes 
and makes one wonder whether a mean-field solution can 
be valid at all in any finite dimension. 

In this paper we revisit the self-consistent mean-field 
solution of Langevin equations with MN and report on 
a new previously overlooked regime. After that, we dis- 
cuss under which circumstances such a solution breaks 
down (Ginzburg criterion |28l l29jl. Also, we present a 
fluctuating-solution aimed to extend the self-consistent 
one by allowing for fluctuations of the average field. Fi- 
nally, we present some speculations on the issue of the 
existence of an upper critical dimension for this type of 
systems, i.e. if there is a finite space dimensionality above 
which the non-fluctuating solution holds or not. 

The paper is organized as follows. In Section [HI 
we outline a self-consistent solution of the MN equa- 
tion defined in an infinite fully connected lattice. Dif- 
ferent regimes are found, corresponding to weak-noise, 
intermediate-noise, and strong-noise respectively. By 
constructing a Ginzburg criterion we will show how the 
strong-noise solution is intrinsically unstable as soon as 
fluctuations are taken into account. In Section II I II we 
study the solution in finite lattices (fluctuating solution) 
and show by means of numerical simulations how the 
mean-field solution breaks down in the vicinity of the 
critical point. Section llVl is devoted to the study of mul- 
tiplicative noise on random regular graphs. This is done 
in order to ascertain the type of behavior in the thermo- 
dynamic limit when the connectivity of each site remains 
finite and, therefore, in order to obtain some insight on 
the behavior of the MN equation in large but finite space 
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dimensions. The final section contains a summary of the 
findings along with some concluding remarks. Finally, a 
general scaling theory is presented in an appendix. 



II. NON FLUCTUATING SOLUTION 

We summarize and extend the mean-field solution ob- 
tained by Birner et al. 0| for the MN equation as de- 
scribed in H,E|: 



4>i = afa - b(f> P+L + DV 2 (j)i + fcrno 



(1) 



where 77 is a Gaussian delta correlated noise, a, b, D, 
and a are constants, and <fi is the local order-parameter 
field, describing the physical density under study. In or- 
der to study its mean-field solution, we consider the case 
of global coupling, i.e. define the system on a fully con- 
nected (or complete) graph 



N 



3^ 



i) = Mi 



(2) 



The associated Fokker-Planck equation for Eq. in 
the Ito sense can then be solved in the stationary case |9j, 
[H[l3, assuming M; = M, 



exp 



2b(f>P 2MB 



pa* 



(3) 



This is a power-law distribution function with two cut- 
offs: an upper one coming from the non-linear saturation 
term and a lower one generated by M which acts as a 
constant external field. 

In order to proceed further M is taken equal to its 
average over realizations, ro, which is determined self- 
consistently by imposing [3, HH 03 



(4>) = 



J d(f>(f>R((p; m) 
J d<j)R(<p; m) 



(4) 



We will denote in the following this solution, valid for 
N = 00, as the Non Fluctuating (NF) solution. Equa- 
tion J3J has always the trivial solution m = 0, stable for 
a < a c with a c = 0. At the critical point, a c , a stable 
solution with m > appears. Introducing the distance 
from the critical point e = a/(a 2 /2) and the reduced 
variable s = 2D /a 2 , we can rewrite the probability dis- 
tribution as 



R{(j)\m) — cj) a exp 



bscj) p ms 



pD 

Introducing the notation 



d(j) (f> k R s (<j)\m) 



a = 2 + s-e. 

(5) 

(6) 



the first moment m — / s ,i// s ,o is easily computed, yield- 
ing the relation between m and e in the active region 



m min[ S -e,p] _ £ 



(7) 
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FIG. 1: Plot of the second, third and fourth moment as a 
function of the first one in the NF strong-noise regime, for a 
system with p = 2, N = 5000 and s — 2/3. All moments scale 
in the same way. 



so that m ~ e^ 1 with (3\ = max[l/s, In this way, 

two different regimes were found by Birner et al.: one is 
universal (weak-noise) and the other one is not (strong- 
noise) with the order-parameter exponent changing con- 
tinuously with s. These two regimes are the analogue 
of the well-known weak-coupling (Gaussian) and strong- 
coupling (non Gaussian) regimes of KPZ dynamics 19]. 

Let us now go beyond the results in |12( by computing 
higher moments = (<fi k ) = I s ^/I s .o- The normaliza- 



tion i, n scales as 



,-l-a+e 



Is,k — I, 



B k m 



-s+k-l+e 



and to leading order 



D k m p ~ s+k - 1+f -. 

While the third term is always sub-leading, it depends 
on k which of the other two dominates. For s > k — 1 + e 
the first one is dominant, so that 



m k 



1 



^k — l — s-\-e 
„-l-s+c 



(9) 



and for s < k 
that 



1 + e the second is the leading one, so 



,1+s- 



(10) 



Hence, the regime with non universal first moment is 
quite rich. For l/p<l/s<l the first k moments exhibit 
standard scaling, as long as k < s + 1 — e, while all others 
scale with the same exponent 1/s+l. We call this regime 
intermediate-noise; it does not have a clear analogue in 
free KPZ-like interfaces. 

When 1/s > 1 full multi-scaling occurs: all moments, 
except the first, scale with the same exponent, 1/s + 1. 
This is the strong-noise regime; numerical evidence of 
this multi-scaling is provided in Fig. ^ and is simi- 
lar to the known phenomenology of single-site (zero- 
dimensional) MN equations [IJlI- 

^From naive power counting (usually expected to re- 
produce mean-field exponents) performed on the MN 
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e 1 1 / s , which diverges 
e. in any case cite |31|. 



equation, the naive time scale T should scale as T _1 ~ 
oae and, therefore, we expect the critical time-decay ex- 
ponent of the order parameter, 9\ (or simply 9) defined by 
((</>(£, e = 0)) ~ t~ dl ) to coincide with (3i in mean-field. 
The same property remains valid for higher moments, so 
the multi-scaling property of the static exponent is trans- 
lated into multi-scaling of the decay exponents [30| . 

Let us now construct a Ginzburg criterion within the 
fully connected lattice, to see under which circumstances 
the previous approximation ceases to be sound. To do so 
we just need to compute the ratio rr ^i 2 [28l l29| . When- 
ever r diverges fluctuations are expected to play a sig- 
nificant role, breaking down the self-consistent solution 
which ignores them, as soon as they are taken into con- 
sideration. From the previously reported scaling for the 
strong-noise solution one has r 
at the critical point if 1/s > 1, 

Therefore, the previously reported strong-noise mean- 
field solution is fully valid only in the strict N = oo 
limit, when fluctuations can be safely discarded owing 
to the law of large numbers. In general, it is expected 
that such a strong-noise solution will break down as soon 
as fluctuations are somehow taken into account. If, for 
example, the MN equation is defined on the top of a 
d-dimensional lattice (every site sees a finite number of 
others) the NF strong-noise solution is not expected a 
priori to be valid. We will discuss this aspect more in 
detail in the forthcoming sections. On the other hand, 
for the weak- and the intermediate-noise regimes, we have 
Gaussian scaling of the lowest moments, and therefore r 
converges to a constant at the critical point. 



III. FLUCTUATING SOLUTION 

In the previously reported approach the crucial step 
is the replacement of M (the mean value seen by any 
arbitrary single site) by a fixed non-fluctuating m. As 
said before, owing to the law of large numbers this is ex- 
act if N = oo whatever the probability distribution of 
the neighboring sites, but one can wonder whether this 
replacement is acceptable if a finite system size A is con- 
sidered. In other words: what is the critical behavior for 
large but finite values of iV? This is a perfectly sensible 
question, since for MN a sharp phase transition is well 
defined for any value of A, even for N = 1 

In this section we relax the condition of a fixed value 
for M, by allowing the spatial average value of the 4> 
field in the stationary state to fluctuate in time, sam- 
pling some probability distribution Q(M) to be deter- 
mined self-consistently. 

Assuming that M changes in a time scale much larger 
than the characteristic time scale for <j>, we can still solve 
the Fokker-Planck for a fixed value of M, and therefore 
Equation © still holds, but now it has to be interpreted 
as a conditional probability R(<fi | M). The full distribu- 



tion P(4>) is given by the convolution: 

/>oc 

P(<P) = / dMR{4>\M)Q{M). 



(11) 



In this case, the self-consistent Equation J3J has to 
be replaced by a self-consistent functional equation for 
Q(M). 



dM M Q{M) 



J ty <t> J dMQ{M)R{4> | M) 
J d(j> J dMQ(M)R(cj> | M) 



(12) 



By solving this functional equation one could obtain the 
full solution as in the NF case, within the slow changing 
M approximation. 

Instead of solving numerically such a self-consistent 
equation we now try to write down an evolution equa- 
tion for M(t), from which Q(M) follows. For that, we 
sum Eq. Q over i and divide it by A, giving 



M = aM-b 



Vi 



(13) 



This turns out to be an effective MN equation for the field 
M in zero dimensions. Let us justify this statement. The 
second term on the r. h. s. is perfectly analogous to the 
nonlinear term in Eq. ft): it just introduces an upper 
cutoff M u in the distribution of the field M. The last 
term on the r. h. s. is less trivial: it can be rewritten 
as o'j^^2 i ^r) i M, hence in the form of a multiplicative 
noise rj'M. The average value of rf is clearly zero, since 
77 has random sign. Its second moment, that we denote 
as cr 2 (A), is 



M 



Vi 



(14) 



The variables r\i and <f>i are uncorrelated. The r\i are 
normally distributed. In practice, to compute cr 2 (A) 
we must evaluate the sum of A variables <f>i/M, dis- 
tributed according to R(<f>/M\M) with random signs. 
The distribution R(<f)/M\M) is a power-law with expo- 
nent a = 2 + s — e, lower cutoff in 1 and upper cutoff in 
1/M. 

We now discuss the properties of this solution depend- 
ing on whether we are working in the strong-noise regime 
or not. 



A. Weak- and Intermediate-Noise 

Let us consider first what happens when 1/s < 1, i.e. 
in the weak and the intermediate noise regimes. In these 
cases, the exponent a of the distribution R is larger than 
3, for a sufficiently small. Adding A such variables with 
random signs is equivalent to adding Gaussian variables 
HI. Then 



A. 



(15) 
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so that 



a 2 (N) ~ o- 2 /7V. 



(16) 



We can conclude that the dynamics of M is governed by 
an effective equation with MN in zero dimensions with 
a noise rj with renormalized variance a 2 /N . ft is clear 
that 77' is correlated in time, but its finite correlation-time 
can be eliminated by suitably rescaling the time variable. 
The distribution for M is then given by the solution of 
the zero-dimensional MN Q 



Q(M) = M 



-2[l-a/cr 2 (AT)] 



exp(-M/M„) 5 



(17) 



As expected, this distribution function is a solution of 
the previously written functional self-consistent Equa- 
tion jnji. 

For finite N the system undergoes an absorbing phase 
transition for a finite value of the control parameter 



a c (N) 



a(N)< 



a 

2N' 



(18) 



The exponents for the transition with finite iV are given 
by the values for zero-dimensional MN transition 0, Q : 
ft = 1, and 9 = 1/2. 

Notice that this transition is qualitatively different 
from that occurring in the thermodynamic limit. Here 
the transition takes place because the distribution of 
M becomes non-normalizable due to the divergence for 
M — > 0. In the thermodynamic limit, instead, the dis- 
tribution Q(M) remains normalizable (Gaussian) at the 
transition and criticality comes from the peak position 
M u moving toward zero (see the appendix, where a co- 
herent general scaling picture is presented). 

These two regimes are therefore distinguished by the 
presence of a (Gaussian) peak for finite M (the NF limit) 
or a broad distribution with a power-law divergence for 
M — * 0. The crossover occurs where the power-law ex- 
ponent in Eq. I|17l) changes sign, i.e. for 



5(JV) = a 2 (N) = 2a c (N). 



(19) 



The behavior for fixed N can then be summarized as 
follows. For a ^> a(N) the system exhibits the criti- 
cal behavior of the weak noise N — 00 solution. For 
< a — a c (N) <C a c (N) it behaves as it was zero- 
dimensional. In the latter case, the crossover can also 
be observed with fixed a by looking at the temporal evo- 
lution of the first moment m: at short times the NF 
solution is followed, crossing-over later to the asymptotic 
zero-dimensional scaling. 

We have checked the correctness of this scenario by 
means of numerical simulations of Eq. Q with p = 2. 
From the numerical point of view, the first problem is the 
determination of the critical point a c (N). The criterion 
we have chosen is based on the way the first moment 
m(t) decays in time. a c (N) is the value separating a 
concave behavior (for a > a c (N)) from a convex one (for 
a < a c (N)). In this way, we obtain the values plotted in 




FIG. 2: Plot of a c (N) versus N for three values of 1/s. Thick 
lines are power-laws fits with exponents —1 for 1/s = 1/8, 
—0.65 for 1/s — 3/2 and —0.57 for 1/s = 2 (to be compared 
with the output of Equations 1181 and 12111 ) . 



Fig. The behavior found in the weak-noise regime is 
in perfect agreement with the predicted 1/N behavior. 

With p = 2, the temporal behavior in the weak- noise 
case is the same (9 = 1/2) both in the thermodynami- 
cal limit and for the effective zero-dimensional behavior 
valid at finite N, not allowing to distinguish between the 
two regimes. In Fig. [3] we analyze the behavior of the 
stationary value of the first moment m. In the main 
part of the figure we observe that, for N = 5000, the 
system follows very accurately (for the values of a con- 
sidered) the decay with ft = 1/2, expected in the ther- 
modynamic limit. For N = 10, instead, a singularity 
is observed for finite a. If we plot m as a function of 
a — a c (N) (Fig. inset) both behaviors can be observed: 
for a - a c (N) > a c (N) w 0.014 the scaling is of the NF 
type (ft = 1/2). For a - a c (N) < a c (N) we observe the 
zero-dimensional exponent ft = 1. 

A direct validation of Eq. (|17fl is provided by Fig. 0] 
where the distribution of the self-consistent field M ex- 
hibits a peak in M u (a) for a — a c (N) ^> a c (N) while a 
power-law divergence at zero develops for a — a c (N) <C 
a c (N). 



B. Strong-Noise 

Let us consider now the strong- noise case 1/s > 1. The 
general picture is similar to the one previously described, 
with the main difference that the exponent a in the dis- 
tribution R(4>/M\M) is now smaller than 3. From the 
theory of Levy-stable distributions [32] we know that the 
sum of N variables distributed as a power-law with ex- 
ponent 1 + [i with 1 < p, < 2 and an upper cutoff 1/M, 
scales as N 1 ^ for N <C M ~ p while it behaves in a Gaus- 
sian way, iV 1 / 2 , for N ^> M - ^. This implies that, for 
fixed TV, t)' is a power-law distributed noise with expo- 
nent 1 + fi for M <C N^ 1 ^, and a Gaussian noise with 
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FIG. 3: Main: plot of m vs a for 1/s = 1/8 (weak- noise). 
Inset: the same data of the main part with N = 10, plotted vs 
a — a c (N) and showing the NF behavior (m ~ (a — a c (A r )) 1 ^ 2 ) 
and the zero-dimensional one (m ~ (a — a c (N))). Solid lines 
are guides to the eye. 



FIG. 5: Main: plot of m vs a for 1/s = 2 (strong-noise). 
Inset: the same data of the main part with N = 10, plotted vs 
a — a c (N) and in agreement with the zero-dimensional linear 
behavior. Solid lines are guides to the eye. 




FIG. 4: Plot of Q(M) vs a for 1/s = 1/8 (weak noise), dis- 
playing the crossover between a Gaussian around M u > and 
a power-law divergence in 0. The value of the critical point is 
a c (N — 10) ~ 0.014. The values of a are, from right to left: 
0.4, 0.1, 0.05, 0.03, 0.02, 0.155. 



a 2 (N) ~ a 2 /N for M > iV 
fi = a — 1 = 1 + s — 2a/a 2 . 



-l/ix 



In the present case 



Eq. Ijl^jl describes now a zero-dimensional MN with a 
rather exotic noise, whose distribution depends on M. 
We have no clear theoretical understanding of such a 
model. In principle it could give rise to completely new 
and non-trivial critical features. However, as shown nu- 
merically below, it turns out to behave asymptotically 
as the zero-dimensional case with standard noise. The 
only change is in the position of the critical point a c (N). 
Accordingly, we assume, rather crudely, that the power- 
law noise does not change the behavior of the zero- 



dimensional MN, except for the form of <r 2 (N) 



<J 2 (N) ~ cr 2 ^ 1 /^-!)-!) = er^Il/U+s^a/CT 2 )-!] 

(20) 

The distribution of M is then given again by Eq. I|17fl . 
with the additional complication that a 2 (N) depends on 
a. The critical point is determined by the implicit con- 
dition 



a c {N) 



a 2 (N) 



a 2 N 



1 



(1 



2a c (N)/a 2 ) 



- 1 



(21) 

If, as a first approximation, we neglect the dependence 
of a on a, i.e. we take a — 2 + s, we obtain a c (N) ~ 
jy-2s/(i+s) ^ n0 |. £ ar f rom t ne numerical results of Fig. [21 

Again for a > a(N) = 2a c (N) the exponent of the dis- 
tribution Q(M) becomes negative and there is a crossover 
to the NF limit. In Fig. [S] we plot m versus a in the 
strong-noise regime. Again the value /3i = 1/s, valid in 
the thermodynamic limit, is observed for large N, as the 
critical point and the crossover are very close to 0. Near 
the transition the exponent is the zero-dimensional one, 
Pi = 1, indicating that the non trivial noise does not 
modify the zero-dimensional critical behavior. 

As in the other case, it is interesting to look also at 
the temporal evolution of m. In the strong-noise regime, 
since the exponent 8 is different in the NF solution and 
in the zero-dimensional regime, the crossover between 
the two regimes results in a crossover in the decay of 
m(t). This is evident from Fig. The effective exponent 
switches from a value close to the prediction 9 = 1/s = 2 
for short times, to the zero-dimensional value 9 = 1/2, 
for longer times. 
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FIG. 6: Plot of m vs t for 1/s = 2 (strong- noise) and a ps 
a c (iV) for (top to bottom) N = 50, N = 500, and N = 50000. 



C. Discussion 

Let us underline the non-commutativity of the limits 
a — > a c and AT — > oo. It is only when the thermody- 
namic limit is taken first and a homogeneous mean-field 
M is considered, that the NF solution is recovered. On 
the other hand, if one takes first the limit a — > a c {N) 
for a generic finite value of iV, the zero-dimensional so- 
lution always dominates the scaling nearby the critical 
point, no matters how large N . Therefore, the inclusion 
of fluctuations has a dramatic effect on the NF solution. 

Such a conclusion holds for all the regimes considered: 
weak, intermediate, and strong-noise. This is due to the 
fact that the two aforementioned limits do not commute, 
implying the presence of a non-analyticity of the most 
general solution in a neighborhood of the critical point. 
This is analogous to the observation of Gaussian scaling 
in standard phase transitions whenever the system size 
is not infinite; it is only when the thermodynamic limit 
is taken that the true asymptotic scaling emerges. In 
the case studied here the role of the Gaussian scaling 
is replaced by a MN zero-dimensional non-trivial scaling. 
This breaking down of the thermodynamic- limit behavior 
for finite systems has, in principle, nothing to do with the 
previously constructed Ginzburg criterion which leads to 
a breakdown of the mean-field solution only at the strong- 
noise regime. 

In the appendix we present a general scaling theory 
accounting in a compact form for all the previously dis- 
cussed phenomenology. 




FIG. 7: Plot of a c (N) versus N for 1/s = 2 (strong-noise) for 
the Fully Connected and the RRG. 

any given one is finite (2d for a hyper-cubic lattice) we 
should first answer the following question: does the NF 
behavior emerge because the size of the system goes to 
infinity or because the number of nearest neighbors di- 
verges? (Note that in the fully connected graph these 
two limits coincide). In order to clarify this point we 
have considered the MN on a connected Regular Random 
Graph (RRG), where each site has fixed degree k > 2 and 
random connections (for k = 2 we have a one-dimensional 
lattice). In this the number of neighbors is finite 

for each site, we expect fluctuations in M and therefore 
a possible breakdown of the NF solution (at least in the 
strong-noise limit), even if the large system size limit is 
taken. 

However, numerical results disprove such a conjecture, 
as we show in what follows. We have performed simula- 
tions of a system with k = 10 and growing N. It turns 
out that the position of the critical point a c (N) depends 
on N and does not reach a finite value dependent only 
on k (see Fig.[7J). Its behavior is not very different from 
what occurs on the fully connected system (see Eq. 1211 . 

Moreover, one can monitor the temporal behavior of 
m(t) for a w a c (N) (Fig.|8J. One finds a crossover from a 
NF behavior at short time to a zero-dimensional behavior 
at longer time, exactly as for the fully connected graph. 

We conclude that the thermodynamical limit is de- 
scribed, also for a regular random graph with finite con- 
nectivity, by the NF behavior. The zero-dimensional be- 
havior holds only as long as criticality is studied for finite 
system sizes. A discussion of these facts is presented in 
the next section. 



IV. FINITE CONNECTIVITY: RANDOM 
REGULAR GRAPHS 

In order to shed some light on the question of whether 
the NF solution holds or not for an arbitrarily large space 
dimensionality, d, in which the number of sites seen by 



V. SUMMARY AND CONCLUSIONS 

In this paper we have investigated the properties of 
the mean-field solution for systems with multiplicative 
noise. First, we have revisited the self-consistent solu- 
tion obtained in Ref. 01 by assuming that M is a ho- 
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FIG. 8: Plot of m vs t for 1/s = 2 (strong- noise) and a ~ 
a c (N) for (top to bottom) N = 50, N = 500, AT = 5000 and 
N = 50000 on the RRG. 



mogeneous non- fluctuating field. We have shown that 
three regimes can actually be identified depending on the 
noise amplitude 1/s = a 2 / (2D). In the weak-noise case 
(1/s < 1/p) all moments of the order parameter distribu- 
tion obey ordinary Gaussian scaling. In the intermediate 
regime, (1/p < 1/s < 1), Gaussian scaling holds only up 
to a certain order, with multi-scaling emerging for higher 
moments. In the strong-noise regime (1/s > 1) all mo- 
ments scale with exponent l/s + 1 except for the first one 
that goes as e 1 ^. 

Complete absence of fluctuations for the spatial aver- 
age M holds only in the thermodynamic limit, with a 
strictly infinite number of nearest neighbors. In this pa- 
per we have relaxed the condition that M can assume 
only one (self-consistently fixed) value and let it fluctu- 
ate: this is equivalent to considering a fully connected 
graph with a finite number N of nodes. 

When N is finite M fluctuates and the critical behav- 
ior is not described anymore by the NF thermodynamic 
limit, in none of the three previously described regimes. 
While in the NF case, criticality is governed by fluctua- 
tions of the 4> field in a single realization around its mean 
value, for finite N what actually matters are the fluctu- 
ations of the self-consistent field M. They are described 
by an effective equation for a zero-dimensional system 
with multiplicative noise, which exhibits a critical behav- 
ior distinct from the one displayed in the thermodynamic 
limit. 

The crossover between the two types of behavior is well 
described by the zero-dimensional description: both the 
crossover point and the critical point location depend on 
N so that for any finite value of N there is a finite interval 
of values of the control parameter where the effective 
zero-dimensional behavior is observed. 

Finally, by means of numerical simulations on con- 
nected Random Regular Graphs, we have also shown 
that all these results still hold (at least qualitatively) for 
random regular graphs with finite connectivity. This is 



somehow counterintuitive, and seems to contradict the 
previously introduced Ginzburg criterion for the strong- 
noise regime. Indeed, given that the fluctuations of M 
are predicted by the Ginzburg criterion to be relevant 
in the strong-noise regime of the NF solution when the 
number of nearest neighbors is not infinity, it is hard 
to understand, why the solution in the random regular 
graph, in which Mi is the fluctuating average taken over 
the k neighbors of any given site i, behaves so similarly 
(for any value of N) to the solution in the fully connected 
network in which M = Mi is fluctuation-less. 

We do not have a clear explanation of this fact, but 
we believe that this is so because of the small- world |23] 
nature of the RRG topology. The small-world property 
implies that one can reach any arbitrary node starting 
from a generic site with a small number of steps following 
network links. In this way, every site is nearby any other 
one, making it difficult to create "local patches" with 
an over-density or sub-density, which would give rise to 
inhomogeneities and a broad field distribution. At this 
point, it would be interesting to study the MN equation 
on a Cayley tree to see if, by introducing well defined 
spatial neighborhoods, the previous results and interpre- 
tation are sustained. We are presently analyzing such a 
problem, but prefer to leave the delicate issues involved 
in such a study for a future publication. 

Many interesting questions remain to be answered. 
The main one is whether the obtained NF solution in the 
strong-noise regime is valid in arbitrarily large but finite 
physical dimensions, or whether it emerges only in d = oo 
(which is obviously related to the existence of a finite up- 
per critical dimension for MN and KPZ problems). The 
Ginzburg criterion for the strong-noise solution seems to 
point out to the second possibility, while the fact that 
in the finite-connectivity random regular graphs the NF 
solution emerges in the thermodynamic limit, could be 
interpreted as supporting the first one. Further analysis 
along these lines is left for future work. 

It remains also to be understood which is the fate of 
the intermediate-noise regime once the MN equation is 
embedded into a d-dimensional lattice. Does it simply 
disappear? Does it survive, introducing some type of 
anomalous effect for high-order moments? Does it have 
any analogous in KPZ-like systems? 

It is our hope that this work will stimulate future re- 
search in this exciting and ever surprising field of systems 
with multiplicative noise. 
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Appendix: General scaling of the distributions 

It is possible to summarize all the scaling regimes de- 
scribed in the paper in a general scaling form of the 
distributions. This describes the crossover between the 
zero-dimensional case and the thermodynamic limit as 
N grows. The equations obeyed by the distribution 
R((f>\M; a, N) of the <fr field conditional to a value M and 
by the full distribution P(4>; a, N) are 



f R(<f>\M;a) oc cxp ( 
\P(4>;a,N) = J* °° dMR((j)\M;a)Q(M;a,N). 



bs<t> p 
pD 



(22) 



The solution of these equations found assuming no fluc- 
tuations for M is 



' R((f>; M) cx 0-«exp(-^-^) 

Q(M;a,N) » 8[M-m(a)] ( 23 ) 
k P(<f>;a,N) = R(4>\M)\ M=m(a) 



A more general explicit ansatz for the solution to this set 
of equations, valid for generic N, is 



R{<t>\M) oc M"- 1 
Q(M; a, N) « M 



P(cj>;a,N) 



*-r <p< M u (a) 



-bs<f,P 

pD 



> M u {a) 



where a' = 2[1 - a/<r 2 (N)]. The form of Q(M) encodes 
the different critical behaviors. In the NF limit (N going 
to oo first) a 1 is negative. In this way Q(M) tends to 
5(M — M u ); moments go to zero because M u goes to zero 
as a 13 . In the opposite nontrivial limit (e going to zero 
with N fixed) a' is positive and grows up to 1 for a — > 
a c (N). In this case, the moments go to zero because the 
distribution Q(M) develops a peak in zero. The critical 
behavior of the moments in this limit is governed by the 
way a' goes to 1. The origin of the two limits (e — > 
and N — > oo) non-commutativity has its roots in the 
non-analytical form of Q(M at M = 0. 



(24) 
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